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ABSTRACT 

We observed the blazar 0716+714 with the VLBA during its active state in 2003-2004. In this paper we discuss multi-frequency 
analysis of the inner jet (first 1 mas) kinematics. The unprecedentedly dense time samphng allows us to trace jet components without 
misidentification and to calculate the component speeds with good accuracy. In the smooth superluminal jet we were able to iden- 
tify and track three components over time moving outwards with relatively high apparent superluminal speeds (8.5-19.4 c), which 
contradicts the hypothesis of a stationary oscillating jet in this source. Component ejections occur at a relatively high rate (once in 
two months), and they are accompanied by mm-continuum outbursts. Superluminal jet components move along wiggling trajectories, 
which is an indication of actual helical motion. Fast proper motion and rapid decay of the components suggest that this source should 
be observed with the VLBI at a rate of at least once in one or two months in order to trace superluminal jet components without 
confusion. 
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OO 1. Introduction stead of travelling along the jet, components are "oscillating" 

around their mean position due to the jet precession. Component 

- The blazar S5 0716+714 is one of the most active sources of ^^^^^^^ calculated taking position angle changes into account, 

. ^ ; Its class. It IS highly variable on time scales from hours to ^^^g^ ^^^^ 5 c to 10c (Britzen et al. 2006). 

^ months at all observed wavelengths from radio to y-rays Intra- ^^.^ ^^.^^ ^ ^^^.^^ ^j^^^ ^^^^^.^^ 

_ day variabihty (IDV) was detected m the source in the op- parsec-scale structure of 0716+714 during the 

P. , tical rmllim etrg^nd radio bands (e g. | Montagm et a . | |20Q6t ^^^^^ -^^ 2003-2004. Here we present the kinematics of 

l ^g'd»^t al.| \2m and references therein). The IDV in the -^^^^-^^ 0716+714. In Paper II we will discuss the suit- 

differe nt bands of electromag n etic radiation is str o ngly corre- ^^-^^ ^-^^^.^^^ ^^^^ ■ ■ ^^^^^^^ f^^. restoration of the 

lated (Ouiii-enbach et al 1991; Wagn er et al. || l996|; I Stalm et al.| ^^^^.^^ ^t^ucture in the case of a very smooth brightness distri- 

L2006, Fuhrmann et al. 2008), which suggests its mtrinsic on- ^^^.^^ ^^^.^^ 3^^^.^^ polarisation, quantitative analysis of 

gm. Redshift ofdie source is 0.31, as recently estimated by ^j^^ evolution and spectral properties will be pub- 

i Nilsson et al l (12008), who observed the host galaxy in the I-band j^^j^^^ subsequent papers. Here we use a standard cosmo- 

durmg the quiescent state of the active nucleus. j^^.^^j ^^^^j ^^^j^ ^ ^^^^^^^ ^^j^^^ p^^^^j^ ^^^^^^^^ 

Very long baseline interferometry studies (VLBI) show a ^^^71 s-'Mpc', energy density of matter Q„,=0.3, and cos- 

core-dominated jet pomtmg to the north (j Jorstad etalj |2001t mological constant energy density Qa=0.7. 
iBach et al.ll2005l), and VLA data shoy y a halo-like jet misaligned 
with it by ~ 90° (I Wagneret ai1ll996h . There are several scenar- 
ios for the milliarcsecond scale jet kinematics. The most recent 
studies propose two very different interpretations for the long 

time series (~10 yr) of centimetre VLBI observations. One see- 0716+714 was observed with the Very Long Baseline Array 

nario consists of fast superluminal compo nents moving along (VLBA). Five observations, separated by approximately a month 

the jet with speeds ranging from 5 c to 16 c ("Jorsta d et al1l200U (2004, Feb 10 (A), May 03 (B), Jun 18 (C), Jul 29 (D) and Aug 

[Bach et al. 2005). An alternative interpretation suggests that in- 29 (E)), cover a time period of the source's active state. This is an 

unprecedentedly high time resolution for this source. Each ob- 

Send offprint requests to: E. Rastorgueva, e-mail: elirasOutu. fi servation lasted for nine hours and utilised five frequencies (86, 



2. Observations and data analysis 
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Fig. 1. Contour maps of 0716+714 at 43 GHz for five epochs convolved with the same beam. Contour levels are 7, 14, 28, 56, 1 12, 
224, 448, 896, 1792 mJy/beam. The distance between maps in horizontal direction is proportional to the time separation between 
epochs. Grey circles represent Gaussian model components. 



43, 22, 5 and 1.6 GHz) in dual polarisation mode. Data were 
correlated at the Socorro VLBA correlator. A priori amplitude 
and phase calibration was performed at Tuorla Observatory us- 
ing standard procedures of the AlPS package. Imaging was done 
using Caltech Difmap package employing CLEAN algorithm. 

2. 1. Model fitting and error analysis 

After the final self-calibrated data and CLEAN image were ob- 
tained, we fitted a source model consisting of Gaussian com- 
ponents to the data in the (m, v)-plane using Difmap procedure 
Modelfit. We used only circular Gaussian components in order 
to reduce the number of free parameters, and we tried to use as 
few components as possible. For all models we took the bright- 
est component at the end of the jet as a reference (the "core") 
and the component positions were measured relative to it. 

The errors in position and flux density of the components 
were calculated using Difwrap package (Lovell 2000), which 
employs Modelfit in a loop for finding maximum possible de- 
viation of the component paramete r from the best-fi t value . We 
used it in a procedure proposed bv ISavolainen et al] (l2006l) . ap- 
plying our own criteria for limiting error determination, which 
will be described in detail in the next paper. 

Expected values of ejection epochs and corresponding errors 
were determined using Monte Carlo simulation, as well as ex- 
pected values of component positions and errors in Ra-Dec co- 
ordinates. Each component was replaced with a Gaussian distri- 
bution of values, representing the errors in that point. A weighted 
least-squares fit was calculated for each set in the case of ejection 
epoch determination, producing a (non-Gaussian) distribution of 
zero crossings. The expected value and standard deviation were 
calculated from these distributions using standard equations. 

2.2. Cross-identification of the model components 

Model fitting of the multi-frequency VLBl data set simultane- 
ously is not a trivial task, because resolution, values of model 
parameter errors, and position of the core all depend on the fre- 
quency. Therefore, we have to treat different bands separately. In 
this paper we discuss the kinematics of the inner jet mainly using 
data obtained at 43 GHz. In addition to that, the data set provides 
models of the same region at 86 and 22 GHz. Even though the 
angular resolution at 22 GHz is only slightly lower (s;0.25 mas) 
than at 43 GHz (^0.15 mas), this difference is important for the 
small inner jet region with fast moving components (see Sect.O. 
Slightly higher resolution is provided by 86 GHz, but the data 
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Fig. 2. Radial distances from the core as a function of time with 
linear fits. Stationary component C8 is plotted with triangles. 
Shaded area represents a Icr error of component ejection time, 
vertical dashed lines mark the beginning of the mm flare. 



quality at this frequency allows us to obtain 86 GHz images 
only for two epochs of our experiment. Therefore, our kine- 
matic model is based on 43 and 22 GHz data (with emphasis 
on 43 GHz), while 86 GHz data is used for a consistency check 
of the models. All 22, 43 and 86 GHz models will be used for 
the subsequent component spectra analysis. 

Models provided by the Modelfit procedure are not a unique 
representation of the data, they are only one of all the possible 
ways to parametrise the image. Therefore, we had to work out 
tests that help us to estimate the goodness of the model fit. We 
consider a model acceptable when the following requirements 
are fulfilled: a) model provides a good fit to the visibility data; 
b) if convolved with the restoring beam, it provides a brightness 
distribution similar to that of the final CLEAN map; c) all models 
at the same frequency are consistent between epochs; i.e., most 
of the components are traceable over time; and d) models at the 
same epoch are consistent between frequencies. 

Condition a) is fulfilled if the model-fitting process was 
carried out correctly and verified by visual inspection of the 
correspondence between the observed and model (m, y)-data. 
Condition b) is checked by simple visual inspection of the re- 
sulting images. Condition c) is verified by comparison of com- 
ponent positions at different epochs and flux density evolution, 
no unphysically rapid changes should be seen. In the presence 
of outward motion, goodness of the linear fit is also a good in- 
dicator For condition d) we use the following rule: components 
from two close frequencies with comparable resolution should 
coincide within positional errors and have flux densities that do 
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Table 1. Proper motions and apparent speeds of components. 



Comp. 


yu, mas/yr 




Ejection 


name 




(z=0.31) 


epoch, yr 


C5 


1 n 1 I n nn 

i.ui±u.uy 


iy.4+i. / 


2003. / /+0.05 




1.11+0.13 


21.3+2.5 


2003.86+0.05 


C6« 


0.88+0.03 (1.04+0.07)* 


16.9+1.6 


2003.99+0.01 




0.98+0.04 


18.9+0.8 


2004.02+0.01 


C7« 


0.44+0.04 (0.44+0.16)* 


8.5+0.8 


2004.08+0.04 


C722 


0.76+0.14 


14.6+2.6 


2004.18+0.04 




0.00+0.01 








0.12+0.02 







* In parenthesis: fit performed using travelled distance 

" Values obtained for 43 and 22 GHz are marked with corresponding 



number 



not require unphysically steep spectral indices (or at least have 
similar flux density-time dependence). Also, if two components 
are resolved at higher frequency, while at the lower frequency 
they appear as one component, the distance between those com- 
ponents should not be greater than the resolution at the lower 
frequency, and their total flux density should be close to the 
flux density of the unresolved one. Our final model set seems 
to meet all those requirements with good accuracy. Therefore, 
our 43 GHz kinematic model is self-consistent over time and is 
also in good agreement with both neighbouring frequencies. 

Cross-identification of the model components between the 
epochs is based on their distance from the core, position on the 
(Ra,Dec)-plane and their flux densities. In this paper, we rely on 
the accuracy of the a priori amplitude calibration of the VLBA, 
which has been p reviously shown to achieve ^5% at 43 GHz 
jSavoIainen et alj |2008). Distribution of the polarised intensity 
along the jet is smooth and uniform, therefore it was not possible 
to use polarisation information to cross-identify the components. 
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Fig. 3. Trajectories of motion of the components identified at 
43 GHz. 



travelled distance. Errors in the radial distance fit parameters 
were calculated using components' positional errors, while trav- 
elled distance fit errors were only calculated from the scatter. 
Therefore the former are substantially smaller. It is clear that 
values of proper motions and apparent speeds obtained by both 
methods agree within eiTors, which indicates that the wiggling 
of the component trajectories (see below) does not influence the 
determination of the component speeds. 

The jet components have wiggling trajectories along the jet, 
with the deviation from the overall jet direction greater than po- 
sitional errors (see Fig.|3j. We have to notice that this wiggling is 
less pronounced at lower frequencies (22 GHz), which suggests 
that the inner jet is resolved in transverse direction at 43 GHz. 



3. Results 

The VLBA mapsof0716+714at43 GHz are presented in Fig . [T] 
The date of each observation is indicated in the plot. Well- 
collimated one-sided inner jet of 0716-H714 spans up to 1 mas 
pointing to the position angle of ~27° and has a smooth inten- 
sity distribution. A total of four components were detected dur- 
ing the observations, and one of them (C8) is stationary. Their 
radial distances as a function of time with error bars and linear 
fits for proper motions are presented in Fig. ID The plotted linear 
fits were performed for 43 GHz data (filled triangles/circles). It is 
clear that most 22 GHz components (open triangles/circles) co- 
incide with those at 43 GHz within one or two sigma. At epochs 
A and B, components C6 and C7, respectively, are not resolved 
from C8. Therefore we took positions of C8 as their positions 
at these epochs. At epoch C, components C7 and C8 are not 
resolved at 22 GHz. Only one moving component (C5) was ob- 
served at all five epochs. At epoch E its flux density faded, and 
it split into two components of comparable flux densities. Their 
positions were averaged in order to perform the linear fit. Based 
on C5, we estimated the component lifetime (from ejection to 
fading away) to be ten months. 

Proper motions, superluminal speeds, and ejection epochs 
for all four components for 22 and 43 GHz are presented in 
Table [T] Components C6 and C7 have substantially curved tra- 
jectories, therefore, for them we present proper motions and 
apparent speeds calculated using linear fit for both radial and 



4. Discussion 

4.1. Jet parameters 

Although kinematics results obtained for 22 and 43 GHz are 
similar, all further calculations were performed for the appar- 
ent speeds obtained for 43 GHz, which seem to us more reliable 
due to the better resolution of images. The estimated apparent 
speeds of the jet components are superluminal with a maximum 
value of 19.4 c. The viewing angle that maximises the apparent 
speed is ^sl = ^rctan(y6app) = 3.0°, coiTesponding to a mini- 
mum Lorentz factor of ysL - 7 mm = (1 + ySnpp)'^^ - 19.5 and 
to a Doppler factor of 5sl ~ 19.4. If 9 differs from 0sl, these 
values could be even higher; however, 19.4 is high enough to ex- 
plain the brightness temperature, estimated during the outbursts 
in the framework of the standard synchrot ron source model wit h 
a power-law electron energy distribution (lOstorero et al.ll2006h . 
Our results indicate slightly hig her flow speeds than previous 
VLBI kinematic studies by Bac h et all ([2005) and agree with 
results of Jorstad et al. (2001). However, radio variability gives 
5vai— 10.9 (Hovatta et al..,200 8). This value could be achieved for 
the observed values of the apparent speed if the viewing angle 
flvar = arctan(2ySapp/(y3app + 5l.„-l)-A.5° , which is lower than 
the maximum value of 5.9°, calculated according to the rela- 
tion sin^njax - 2y6app/(l + Therefore, our results are also 
consistent with calculations based on the source variability. 
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4.2. Correspondence of component ejections to mm flares 

According to lSavolainen et alj (|2002|) . an ejection of superlumi- 
nal components in blazars in most cases coincide with the be- 
ginning of the mm-continuum flare. In September-October 2003 
0716+714 underwent an outburst in mm and radio bands. The 
source had a large flare with several peaks and a complicated 
structure in the millimetre band (T. Hovatta and A. Lahteenmaki, 
priv. comm.). In the period under consideration, four flares were 
registered at 37 GHz, starting at 2003.79, 2003.85, 2003.99, and 
2004.05. As a flare start we took one e-folding time before the 
maximum. Three of them coincide with the ejection times of 
components C5, C6, and C7 (see Fig.|2]and Table[T]i. 

4.3. Wiggling trajectories 

We noticed that trajectories of the jet components are not 
stra ight. The standard shock-in-jet model of AGN (MarscheJ 
[1996 ) suggests straight outward movement of the emitting com- 
ponent along the jet, but several alternative models state that 
emitting plasma clouds may move along helical paths. Although 
we do not have enough data to distinguish between those mod- 
els, it might be interesting to apply them to this source in the 
future. 

A two-fluid model of flie jet (ISol et al.1 119891; 
iDespringre & Fraix-Burnej 1 19971) presents the jet itself to 
be non- or mildly-relativistic, consisting of electrons and 
protons, and not radiating. Emission comes from the relativistic 
"clouds" of electron-positron plasma moving along the helical 
path wrapped around the jet and driven by helical magnetic 
field. Brightness distribution on the VLBI map, dependencies of 
the component parameters on the distance from the core, as well 
as synchrotron spectra for certain geometries, were presented in 
IDespringre & Fraix-Bur"ne3(ll997l) . 

The helical jet model (e.g., iGomez et al.l Il994t 
IVillata & Ferraril 1 19941) presents the jet as an inhomoge- 
neous flow of plasma along the helical magnetic field filaments, 
interacting with the ambient matter. It emits constantly, and flux 
density variability is caused by the change in jet orientation. This 
model, if applied to a number of sources incl uding 0716+714, 
explains its spectral variability dOstorero et al.|[2001i) . 

4.4. Long-term variation in the component speeds and 
ejection rate 

iBach et aTl (|2005[) find that the proper motions of the jet compo- 
nents are almost monotonically slowing down during the period 
fror n late 1986 until 1998 (0.86 to 0.29 mas/yr, respectively), 
and lNesci et al.l (l2005b explain this fact by a slow jet precession. 
Their model predicts a further increase in the apparent speed, to- 
gether with decrease in the source's mean brightness. Although 
our data show w ith a certain ty that the speeds in 2004 are faster 
than reported by Bach et alj (120051) . and in the period from 2004 
to 2008, the source is even slightly brighter in optical than in 
1995 - 2003. It indicates that jet precession is not the only mech- 
anism responsible for the changes in the apparent speed. Also 
the ejection interval has reduced dramatically: if each Gaussian 
component in our model corresponds to a distinct physical fea- 
ture in the jet, the ejection rate after the outburst in 2004 is one 
in 0. 17 years, which is much higher than the long-term mean of 
one in 1.2 years. 

The deviations from the mean ejection rate during the pe- 
riod 1986 - 1998 have occurred after small 15 GHz flux density 
enhancements. One explanation for this behaviour could be that 



the mechanism that normally produces a steady rate of ejections 
is disturbed during an outburst and settles down gradually after 
that. 

5. Conclusions 

We studied the kinematics of the inner jet of 0716+714 at 
43 GHz. Even though the jet has a very smooth brightness dis- 
tribution, we identified four components: three moving and one 
stationary. Components are well-defined at several frequencies 
and traceable over time. Components move along the jet with 
superluminal speeds ranging from 8.5 c to 19.4c, which contra- 
dicts the model of an oscillating jet. 

Ejections of new superluminal components occur at a rela- 
tively high rate, once in two months, and a component life time 
is approximately ten months. Therefore, future observers have 
to note that for accurate structure evolution studies, 0716+714 
has to be observed at a rate of once in one or two months. 
Superluminal component ejections coincide with outbursts in a 
millimetre continuum. 

We also noticed that components move along wiggling tra- 
jectories in the jet, which might be a signature of their helical 
motion. Careful quantitative tests of models that predict physi- 
cally helical component motion (e.g. two-fluid model and helical 
jet model) should be performed in the future in order to under- 
stand nature of the jet in this object. 
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